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In this paper we review the theoretical issues involved in the direct detection of super- 
symmetric (SUSY) dark matter. After a brief discussion of the allowed SYSY parameter 
space we focus on the determination of the traditional neutralino detection rates, in ex- 
periments which measure the energy of the recoiling nucleus, such as the coherent and 
spin induced rates and the dependence of the rate on the motion of the Earth (modulation 
effect). Then we examine the novel features appearing in directional experiments, which 
detect the recoiling nucleus in a given direction. Next we estimate the branching ratios 
for transitions to accessible excited nuclear levels. Finally we comment on the event rates 
leading to the atom ionization and subsequent detection of the outgoing electrons. 
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1 Introduction 

The combined MAXIMA-1 [1], BOOMERANG [2], DASI [3] and COBE/DMR 
Cosmic Microwave Background (CMB) observations [4] as well as the recent WMAP 
data [5] imply that the Universe is flat [6] and that most of the matter in the Universe 
is Dark, i.e. exotic. Crudely speaking and easy to remember one has: 

Vl b = 0.05, CIcdm = 0.30, fl A = 0.65 

for the baryonic, dark matter and dark energy fractions respectively. 

These observations, however, do not tell us anything about the particle nature 
of dark matter. This can only be accomplished through direct observation. Many 
experiments are currently under way aiming at this goal. 

Supersymmetry naturally provides candidates for the dark matter constituents 
[7], [8]- [11]. In the most favored scenario of supersymmetry the lightest supersym- 
metric particle (LSP) can be simply described as a Majorana fermion, a linear 
combination of the neutral components of the gauginos and higgsinos [7], [8]- [15]. 

In all calculations performed so far, the obtained event rates are quite low and 
perhaps unobservable in the near future. So one has to search for characteristic 
signatures associated with this reaction. Such are the modulation of the event 
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rates with the motion of the Earth (modulation effect) and the correlation of the 
observed rates of directionally sensitive experiments with the motion of the sun 
[16, 17]. Transitions to low energy excited nuclear states have also been considered 
[18]. Quite recently it has been found that the detection of electrons, following the 
collision of the neutralino with the nucleus may offer another option [19] to be 
exploited by the experiments. 

2 The Essential Theoretical Ingredients of Direct Detection. 

The possibility of dark matter detection hinges on the nature of its constituents. 
Here we will assume that such a constituent is the lightest supersymmetric particle 
or LSP. Since this particle is expected to be very massive, m x > 30GeU, and 
extremely non relativistic with average kinetic energy T m 50 keV(m x /100GeV), 
it can be directly detected [7]- [27] mainly via the recoiling of a nucleus (A,Z) in 
elastic scattering. In this paper, however, we will consider alternative possibilities. 

The event rate for all such processes can be computed from the following ingre- 
dients: 

1. An effective Lagrangian at the elementary particle (quark) level obtained in 
the framework of supersymmetry as described , e.g., in Refs [15, 24]. 

2. A well defined procedure for transforming the amplitude obtained using the 
previous effective Lagrangian from the quark to the nucleon level, i.e. a quark 
model for the nucleon. This step is not trivial, since the obtained results 
depend crucially on the content of the nucleon in quarks other than u and d. 

3. Nuclear matrix elements [28] — [30], [18], [19], obtained with as reliable as pos- 
sible many body nuclear wave functions. Fortunately in the most studied 
case of the scalar coupling the situation is quite simple, since then one needs 
only the nuclear form factor. Some progress has also been made in obtaining 
reliable static spin matrix elements and spin response functions [30], [18]. 

The calculation of this cross section has become pretty standard. One starts 
with representative input in the restricted SUSY parameter space as described in 
the literature for the scalar interaction [13, 15] (see also Arnowitt and Dutta [25]). 

It is worth exploiting the contribution of the axial current in the direct neu- 
tralino detection, since, in addition, it may populate excited nuclear states, if they 
happen to be so low in energy that they become accessible to the low energy neu- 
tralinos [18]. Models which can lead to detectable spin cross sections have recently 
been proposed [20] [21] [22] [23]. 

Once the LSP-nucleon cross section is known, the LSP-nucleus cross section can 
be obtained. The differential cross section with respect to the energy transfer Q for 
a given LSP velocity v can be cast in the form 

d*(u,v) = ^- [(S S F 2 ( M ) + S spm ^ 11 (^)] (1) 
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where we have used a dimensionless variable u, proportional to Q, which is found 
convenient for handling the nuclear form factors [27] F(u) , Fn(u), namely u = 
! Qo ~ 40 x yl~ 4 / 3 MeV. fi r is the reduced LSP-nucleus mass and b is (the 
harmonic oscillator) nuclear size parameter. Furthermore 

= "'f yD^Mri ^spin = Mr°i/v° Cspin i Cspin = To <S'(' U ) (2) 

3(1 + if) 2 

CTp^o and CTp x o arc the nuclcon cross-sections associated with the spin and the 
scalar interactions respectively and 

s(u) = [Sn (oyf^p\ + 2^fi (o)o 1 (o)|^M + ^(o)) 2 ] (3) 

The definition of the spin response functions Fy, with i,j = 0, 1 isospin indices, 
can be found elsewhere [30] . 

Some static spin matrix elements [30], [28], [27] for some nuclei of interest are 
given in table 1 



Table 1. The static spin matrix elements for various nuclei. For light nuclei the calculations 
are from Divari et al (see text) . For 127 1 the results are from Ressel and Dean (see text) 
(*) and the Jyvaskyla-Ioannina collaboration (private communication^**). For Pb they 
were obtained previously (see text). 



19p 29gj 23 Na 127 /* 127 j** 207p^+ 



[^i(O)] 2 


2.610 


0.207 


0.477 


3.293 


1.488 


0.305 


2.807 


0.219 


0.346 


1.220 


1.513 


0.231 


n (o)fii(o) 


2.707 


-0.213 


0.406 


2.008 


1.501 


-0.266 


mm 


2.91 


-0.50 


2.22 








l^exp 


2.62 


-0.56 


2.22 








tJ-th(spin) 


0.91 


0.99 


0.57 









3 Rates 

The differential non directional rate can be written as: 

dRundir = . m do{u,v)\v\ (4) 

m x ArriN 

Where p(0) = 0.3GeV/cm 3 is the LSP density in our vicinity, m is the detector 
mass, m x is the LSP mass and da(u,v) was given above. 
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The directional differential rate, in the direction e of the recoiling nucleus, is given 
by: 

dRdir = , m \v\v.e ®(v.e) ^- da(u,v) (5) 

m x ArriN 2n 

The LSP is characterized by a velocity distribution. For a given velocity distri- 
bution f(u'), with respect to the center of the galaxy, One can find the velocity dis- 
tribution in the lab frame f(v, ve) by writing v=v+ve , ve=vq+ 
V\. Vq is the sun's velocity (around the center of the galaxy), which coincides with 
the parameter of the Maxwcllian distribution, and V\ the Earth's velocity (around 
the sun). Thus the above expressions for the rates must be folded with the LSP 
velocity distribution. We will distinguish two possibilities: 

1. The direction of the recoiling nucleus is not observed. 
The non-directional differential rate is now given by: 

dRundir \ _ /dR\ _ p(0) m g-^ 
du l~ \dul ra* Am N du 



where 



,cE x f \v\ t/ ^da(u,v) 3 



(3-)= / -p=f(v,v E ) V d A v (7) 
du' J yfjy) K ' du w 



2. The direction e of the recoiling nucleus is observed. 

In this case the directional differential rate is given as above with: 



,.cE. , f v.e Q(v.e) „. ,da(u,v) 

— o{ — v.e)d v 

2tt [i r bv 

To obtain the total rates one must integrate the two previous expressions over 
the energy transfer from Q m in determined by the detector energy cutoff to Q ma x 
determined by the maximum LSP velocity (escape velocity, put in by hand in the 
Maxwellian distribution), i.e. v esc = 2.84 u , vo — 229 Km/s. 

4 Results 

We will specialize the above results in the following cases: 
4.1 Non directional unmodulated rates 

Ignoring the motion of the Earth the total non directional rate is given by 

R = Rt(a,Q mm ) (9) 
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- p(\J) m , >2 / / o\ r S A 2 spin * i 

R = ^mT v ( m^) } + c -" ] 

where i is the modification of the total rate due to the folding and nuclear structure 
effects, t depends on Q m in, i-e. the energy transfer cutoff imposed by the detector 
and the parameter a introduced above. All SUSY parameters, except the LSP mass, 
have been absorbed in R. 

Via Eq. (9) we can, if we wish, extract the nucleon cross section from the data. 
For most of the allowed parameter space the obtained results are undetectable. As 
it has already been mentioned it is possible to obtain detectable rates in the case 
of the coherent mode. Such results have, e.g. been obtained by Cerdeno et al [31] 
with non universal set of parameters and the Florida group [32] . 




Fig. 1. On top:The quantity R, « event rate for Qmin = 0, associated with the spin 
contribution in the case of the A = 19 system (for the definition of the parameters see 
text), bottom: The event rate, associated with the spin contribution in the case of the 
A — 127 system (for notation see our earlier work [18]). 



For the target 19 F are shown in Fig. 1 (top), while for 127 1 the corresponding 
results are shown in Fig. 1. 

4.2 Modulated Rates 

If the effects of the motion of the Earth around the sun are included, the total non 
directional rate is given by 

R = Rt[{l + h(a,Q min )cosa)} (10) 

with h the modulation amplitude and a is the phase of the Earth, which is zero 
around June 2nd. The modulation amplitude would be an excellent signal in dis- 
criminating against background, but unfortunately it is very small, less than two 
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per cent (see table 2). Furthermore for intermediate and heavy nuclei, it can even 



Table 2. The parameters t, h, k, h m and a m for the isotropic Gaussian velocity distribu- 
tion and Qmin = 0. The results presented are associated with the spin contribution, but 
those for the coherent mode are similar. The results shown are for the light systems. For 
intermediate and heavy nuclei there is a dependence on the LSP mass. +x is radially out 
of the galaxy (Q = 7r/2, $ = 0), +z is in the sun's direction of motion (O = 0) and +y is 
vertical to the plane of the galaxy (O = 7r/2, <E> = 7r/2) so that (x,y,x) is right-handed. 
Q m = 0,1/2,1,3/2 means that the maximum occurs on the 2nd of June, September, 

December and March respectively. 



type t h dir k h m a. 



+z 

dir +(-)x 

+(-)y 

-z 

"all LOO 

all 0.02 



0.0068 0.227 1 

0.080 0.272 3/2(1) 

0.080 0.210 (1) 

0.395 0.060 



change sign for sufficiently heavy LSP. So in our opinion a better signature is 
provided by directional experiments, which measure the direction of the recoiling 
nucleus. 



4.3 Directional Rates. 

Since the sun is moving around the galaxy in a directional experiment, i.e. one 
in which the direction of the recoiling nucleus is observed, one expects a strong 
correlation of the event rate with the motion of the sun. The directional rate can 
be written as: 

Rdir = t -^-R[l + h m cos(a - a m tt)] (11) 
= -^-R t{l + h m cos(a - a m it)] 

Z7T 

where h m is the modulation, a m is the shift in the phase of the Earth a and k/(2tt) 
is the reduction factor of the unmodulated directional rate relative to the non- 
directional one. The parameters n , h m , a m depend on the direction of observation: 
The above parameters are shown in Table 2 parameter tdir for a typical LSP mass 
100 GeV is shown in for the targets A = 19. For heavier targets the depend on the 
LSP mass [16]. 
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4.4 Rates to excited states 

Transitions to excited states are possible only for nuclear systems characterized by 
excited states at sufficiently low energies with quantum numbers, which allow for 
Gamow- Teller transitions. One such system is 127 I, which, fortunately, can serve 
as a target for the recoil experiment. 

This nucleus has a ground state 5/2+ and a first excited state a 7/2+ at 57.6keV. 
As it has already been mentioned it is a popular target for dark matter detection. 
As a result the structure of its ground state has been studied theoretically by a lot 
of groups (for references see [18]). We find 0,% = fi? = fi ^i = 0.164, 0.312 for the 
ground state and the excited state respectively. 

In presenting our results it is advantageous to compute the branching ratio. 
In addition to factoring out most of the uncertainties connected with the SUSY 
parameters and the structure of the nucleon, we expect the ratio of the two spin 
matrix elements to be more reliable than their absolute values. Taking the ratio 
of the static spin matrix elements to be 1.90 and assuming that the spin response 
functions are identical, we calculated the branching ratio , which is exhibited in 
Figs 2 and 3. We notice that the dependence on Qmin is quite mild. From Figs 




Fig. 2. On the left we show he ratio of the rate to the excited state divided by that of the 
ground state for 127 1 assuming that the static spin matrix element of the transition from 
the ground to the excited state is a factor of 1.902 larger than that involving the ground 
state, but the spin response functions are the same. Next to it we show the modulation 
amplitudes for the ground and the excited states respectively. The results were obtained 
for no threshold cut off (Q m in = 0) . 




Fig. 3. The the same as in Fig. 2 for a lower energy cutoff of Qmin = WkeV. 
2 and 3 we notice that the relative modulation is higher when the phase space 
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is restricted by Qmin and E exc at the expense, of course, of the total number of 
counts. 

4.5 Detection of ionization electrons 

The experiments may also attempt to detect the ionization electrons, which maybe 
emitted during the LSP-nuclcus collision [19]. This interesting signature can occur, 
both for the coherent and the spin mode, with a branching ratio between 2% and 
8%. Due to lack of space, however, we are not going to elaborate further on this 
point. 

5 Conclusions 

In the present paper we have discussed the parameters, which describe the event 
rates for direct detection of SUSY dark matter. In the coherent case, only in a small 
segment of the allowed parameter space the rates are above the present experimental 
goals [13, 15, 25], which, of course, may be improved by two or three orders of 
magnitude in the planned experiments [33]- [37]. In the case of the spin contribution 
only in models with large higgsino components of the LSP one can obtain rates, 
which may be detectable, but in this case, except in special models, the bound 
on the relic LSP abundance may be violated. Thus in both cases the expected 
rates are small. Thus one feels compelled to look for characteristic experimental 
signatures for background reduction, such as correlation of the event rates with the 
motion of the Earth (modulation effect) and angular correlation of the directional 
rates with the direction of motion of the sun on top of their seasonal modulation. 
Such experiments are currently under way, like the UKDMC DRIFT PROJECT 
experiment [38], the Micro-TPC Detector of the Kyoto- Tokyo collaboration [39] 
and the TOKYO experiment [40] . 

The relative parameters t and h in the case for light nuclear targets are essen- 
tially independent of the LSP mass, but they depend on the energy cutoff, Q m in- 
For Qmin = they are exhibited in Table 2. They are essentially the same for both 
the coherent and the spin modes. For intermediate and heavy nuclei they depend 
on the LSP mass [16]. 

In the case of the directional rates it is instructive to first summarize our results 
regarding the reduction of the directional rate compared to the usual rate, given by 
n/(2n). The factors k depend, of course, on the angles of observation (see Table 2). 
Second we should emphasize the importance of the modulation of the directional 
rates. In the favored direction the modulation is not very large, but still it is three 
times larger compared to that of the non directional case. In the plane perpendicular 
to the sun's motion the modulation is quite large (see Table 2). 

Coming to transitions to excited states we believe that branching ratios of the 
size obtained here for 127 1 are very encouraging to the experiments aiming at 7 ray 
detection, following the de-excitation of the nucleus. 

Regarding the detection of the emitted electrons in the LSP-nucleus collision 
we find that, even though the distribution peaks at low energies, there remains 
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substantial strength above 0.2 keV , which is the threshold energy of a Micromegas 
detector, like the one recently [41] proposed. 

We thus hope that, in spite of the experimental difficulties, some of the above 
signatures can be exploited by the experimentalists. 
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